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The angular and frequency characteristics of the gravitational radiation emitted in collisions of 
massless particles is studied perturbatively in the context of classical General Relativity for small 
values of the ratio a = 2 rs/b of the Schwarzschild radius over the impact parameter. The particles 
are described with their trajectories, while the contribution of the leading nonlinear terms of the 
gravitational action is also taken into account. The old quantum results are reproduced in the zero 
frequency limit lo -C 1/6. The radiation efficiency e = E ra ,d/2E outside a narrow cone of angle a 
in the forward and backward directions with respect to the initial particle trajectories is given by 
e ~ a 2 and is dominated by radiation with characteristic frequency ui ~ 0(l/rs). 
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1. INTRODUCTION 

The problem of gravitational radiation in particle collisions has a long history and has been studied in a 
variety of approaches and approximations. The interested reader may find a long and rather comprehensive 
list of relevant references in |l|, where the emitted gravitational energy, as well as its angular and frequency 
distributions in ultra-relativistic massive-particle collisions were computed. The condition imposed in !l| that 
the radiation field should be much smaller than the zeroth order flat metric restricted the region of validity of 
our approach to impact parameters b much greater than the inverse mass of the colliding particles, and made 
our conclusions not applicable to the massless case. 

However, the problem of gravitational radiation in mass/ess-particle collisions is worth studying in its own 
right and has attracted the interest of many authors in the past as well as very recently. Apart from its obvious 
relevance in the context of TeV-scale gravity models with large extra dimensions [2|, it is very important in 
relation to the structure of string theory and the issue of black-hole formation in ultra-planckian collisions [If, 
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@], Nevertheless, to the best of our knowledge, complete understanding of all facets of the problem is still 
lacking. The emission of radiation in the form of soft gravitons was computed in [BJ in the context of quantum 
field theory, but in that computation the contribution of the non-linear graviton self-couplings i.e. the stress 
part of the energy-momentum tensor, was argued to be negligible. The result of the quantum computation for 
low-frequency graviton emission was reproduced by a purely classical computation in Q, due entirely to the 
colliding particles and leaving out the contribution of the stress part of the energy-momentum tensor. 

In the pioneering work [7} or its recent generalization to arbitrary dimensions fsj, the special case of collisions 
with vanishing impact parameter was studied, with emphasis on the contribution to the radiation of the stress 
part of the energy momentum tensor, leaving out the part related to the colliding particles themselves. In a more 
recent attempt |9J the metric was computed to second order, but no computation of the radiation characteristics 
was presented, apart from an estimate of the emitted energy based essentially on dimensional analysis. More 
recently, a new approach was put forward for the computation of the characteristics of the emitted radiation 0 , 
based on the Fraunhofer approximation of radiation theory. However, this method cannot be trusted at very 
low frequencies uCl jb and, furthermore, it ignores the non-linear terms of the gravitational action, which are 
expected to be important in the high frequency regime. Thus, we believe it is fair to conclude, that the issue of 
the frequency and angular characteristics as well as the efficiency of gravitational radiation in ultra-relativistic 
particle collisions is not completely settled yet. 

The purpose of the present paper is to extend the method used in :l| to the study of gravitational radiation 
in collisions of massless particles with center-of-mass energy 2 E and impact parameter b. The formal limit 
to — > 0 (or equivalently 7 cm — > oo for the Lorentz factor) of the massive case leads to nonsensical answers for 
the radiation efficiency, i.e. the ratio e = E Ta a/2E ~ {rs/b) 3 7 cm of the radiated to the available energy, the 
characteristic radiation frequency uj ~ j 2 m /b, or the characteristic emission angle 1 !) ~ l/ 7 cm- The whole set-up 
of the computation in the massive case is special to that case and, consequently, does not allow to extract safe 
conclusions related to massless-particle collisions. In particular, the massive case computation was performed 
in the lab frame, the choice of polarization tensors was special to the lab frame, while, being interested in ultra- 
relativistic collisions, we organized the computation of the energy-momentum source in a power series of the 
Lorentz factor 7 . Here, we shall deal directly with massless collisions in the center-of-mass frame and correct the 
above inadequacies of our previous results. We shall study classically the gravitational radiation in the collision 
of massless particles using the same perturbative approach as in jlj. The scattered particles will be described 
by their classical trajectories, eliminating potential ambiguities in the separation of the radiation field from the 
field of the colliding particles, inherent in other approaches. Furthermore, at the level of our approximation we 
shall take into account the contribution of the cubic terms of the gravitational action to the radiation source, 
which will be shown to be essential for the consistency of our approach. Finally, the efficiency e outside a narrow 
cone in the forward and backward directions will be obtained as a function of the only available dimensionless 
quantity a = 2 rs/b = 8GE/bc 4 , formed out of the four parameters G, E , 6 , c, relevant to the problem at hand. 

The rest of this paper is organized as follows: In Section 2 we describe the model, our notation, the equations 
of motion and the perturbative scheme in our approach. This is followed by the computation in Section 3 
of the total radiation amplitude, i.e. the sum of the local and stress part. Section 4 focuses on the study 
of the angular and frequency characteristics of the emitted radiation in the most important regimes of the 
emission-angle—frequency plane. . Furthermore, in a separate subsection we compare the results of this paper 
to previous work and verify that they are compatible in their common regime of validity. Our conclusions are 
summarized in Section 5, while in three Appendices the interested reader may find the details of several steps 
of the computations and the proofs of basic formulae, used in the main text. 


2. NOTATION - EQUATIONS OF MOTION 

The action describing the two massless particles and their gravitational interaction reads 

S = J e ( (T ) 9 ^{ z (< :T ))z' J '(cr)z I '(a)da- J Ry/^g d 4 x , (2.1) 

where e(cr) is the einbein of the trajectory z ,1 {a) in terms of the corresponding affine parameter a, >c 2 = 167 rG 
and the summation is over the two particles. We will be using unprimed and primed symbols to denote quantities 
related to the two particles. 

Variation of the einbeine gives for each particle the constraint 

g»v(z(cr)) z^{a)z v {a) = 0. (2.2) 

Using the a— reparametrization invariance a —> d = d(a), e(a) -A e(d) = e(cr) dcr/dd we can choose e(a) = 
constant. Furthermore, we can use the remaining freedom of a rescalings to set e = e' on the particle trajectories. 
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Finally, we can shift the affine parameters to set a = 0 = a' at the positions of closest approach of the two 
particles. Before the collision the particle positions are at negative cr and a'. They “collide” when they are at 
cr = 0 = a'. 

For identical colliding particles in the center-of-mass frame we can choose a' = a and, consequently, e = e'. 
With the gauge choice e = constant, the two einbeine are finally determined by the condition 

y/s = E + E' = Jt 00 {x)c1 3 x, (2.3) 


from which 


e = Vs /2 = E , 


(2.4) 


with E the energy of each colliding particle. 

Thus, the particles move on null geodesics, while variation of leads to the particle equation of motion: 


^ {9^z v ) = ^g\v,nZ X z'' 


(2.5) 


and similarly for z At zeroth order in the gravitational interaction, the space-time is flat and the particles 
move on straight lines with constant velocities, i.e . 1 

V = »w; °i M = = (1,0,0,1), °z^ = u^ = (1,0,0,-1). 

The particle energy-momentum is defined by T^ = (—2 /y/—g) SS/Sg , i.e. for each particle 


T^(x) =ej 

At zeroth order, in particular, it is given by 

Qr-pfliy _ 


z^z v 6(x - z(a)) 

7^9 

u^u v S(x — z(a)) 


da. 


V^9 


da 


( 2 . 6 ) 


(2.7) 


and is the source of the first correction h^ v of the gravitational field. Given that °X/„ is traceless, the pertur¬ 
bation satisfies for each particle separately the equation 




whose solution in Fourier space is 




2ttk 


q 2 + i 0 q° 


e**W6(qp) Pli p v and h'(q ) = 


27 TX 


q 2 + iOq° 


( 2 . 8 ) 


e^'^gpVX’ (2.9) 


where = eu^ 1 , p'^ = eit' M , while z^{ 0) = (0, 6/2,0,0) and z'^(0) = (0, —6/2,0, 0). Since the particle momenta 
satisfy p 2 = 0 = p' 2 , the consistency conditions = 0 = 6 '/, are also satisfied to this order. 

In coordinate representation they are 

M*) = = -xeu^u u 5(t - z) $(|r - b/2|) 

(2tt ) 6 J q- 

h 7( x ) = -^^r J ^^e^^)e iq [ r+b / 2 ]=->,eu>U ( t + ^$ ( |r + b/ 2 |), ( 2 . 10 ) 


where $ is the 2—dimensional Fourier transform of 1 /q 2 : 

d 2 q 


1 


$ ( r ) = 7 TI 72 / «-«*' =--In- 


(27r) 


27T r 0 


( 2 . 11 ) 


with r = (x , y) and b = ( 6 , 0 ) is the position and impact vector, respectively, in the transverse x — y —plane and 
ro an arbitrary constant with dimensions of length. 


1 The upper left index on a symbol labels its order in our perturbation scheme. 
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Write for the metric g = rj^ + >f(h^ v + h') and substitute in (12.51) to obtain for the first correction of the 
trajectory of the unprimed particle the equation 


%(cr) = -x ( h' x - -h[ 


Xu, fi 


°z xo z\ 


( 2 . 12 ) 


The interaction with the self-field of the particle has been omitted and due to the primed particle is evaluated 
at the location of the unprimed particle on its unperturbed trajectory. 

We substitute (12.91) into ( 12 . 121 ) to obtain 


(27rj J q L 


Integrating it over a, the first-order correction to velocity is given by 

9^2 


[ d A q 5 ^ e -^ b e ~diu)a 

<7 M 1 
ru H 

U ^ — t-- 

I q 2 

L (q u ) J 


C"h 


(2.13) 


(2.14) 


The integration constants C M are chosen C° = 0 = C z and C x = ex 2 &(b )/2 in order to satisfy the initial 
conditions ^(cr = —oo) = 0 . 

Thus, the components of ^(a) are 


= 


ezt 

(2^)3 

1 


dq 1 


d 2 q 


0 " e iqb e _2i<z a = - ex 2 q>(b) 6{a) 


Y z z (a) = —^ e>c 2 4 > (6) 5(a) 

0 x 2 . 


1 i» = - 

^(cr) = 0 . 
Making use of the formulae [13] 

1 


exr f dq° d ~q 


(2tt) 3 


q 2 


e iqb e -2ig°<T 9 * + C-* = ex 2 


(2.15) 


= ^-i7T ( / J (n ~ 1} (g), 


[x + z 0 ] ra x n (n — 1 )! 


T 


(x + * 0 )”_ 


(k) = 27r 


H) n 

(n- 1)! 


[k9(-k)\ 


n —1 


satisfied by the distributions (x + iO) n and their Fourier transform, respectively, we can express 1 z ft (a ) collec¬ 
tively in the following useful form 




: ■ = -^3 J d A q-^e- iqh e-^° [(qu)u'» - q " 


1 


(gu) + *0 


(2.16) 


which vanish for all er < 0. Indeed, the massless particle trajectories should remain undisturbed before the 
collision. 

Finally, we integrate (12.161) and fix the integration constants so that 1 z fJ ’(a ) is regular and satisfies 
V^er < 0) = 0. We end up with 


V» = [ d 4 q^p- e~ iqb e-W* \(qu)u* - q» 

{Zny J q z L 

or, equivalently, in components 


1 


[(qu) + *0] 5 


(2.17) 


1 2 0 (cr) = i ex 2 $(b) 0(a) = — X z z (a) 


L z x (a) = ex 2 <&'(b) ad(a). 


(2.18) 


From these it is straightforward to reproduce the leading order expressions of the two well-known facts about 
the geodesics in an Aichelburg-Sexl metric, namely 

• The time delay at the moment of shock equal 


At = ex 2 4>(6) = 8 GE In ■ 


r o 
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• The refraction caused by the gravitational interaction by an angle 

2 iii/n 8 GE 
a = ex |$ (6)| = —-— 

in the direction of the center of gravity. 

Clearly, similar expressions to the above are obtained for the primed particle trajectory. For the perturbation 
V M (cr), in particular, we have 


= - 


2 ie>r 


e +^ b e ~i(<i u, y 


{qu')u^ 


. [(qu') + iO}‘ 2 


(2.19) 


To summarize: We have obtained the first order corrections h^, v {x) of the gravitational field, sourced by the 
straight zeroth-order trajectories of two colliding massless particles. It is identical with the leading term of the 
Aichelburg-Sexl metric describing the free particles and it can be shown to coincide with the limit m —> 0 of 
the corresponding field due to massive particles. The perturbations l z^(o) and V^er) of the trajectories of 
the colliding particles in the center-of-mass frame and with impact parameter b were also computed. Finally, 
the known expressions fT^ I for the time delay At and the leading order in rs/b <C 1 scattering angle a were 
reproduced. 

As will be shown in the next section, the arbitrary scale r o in the expressions for h^ and h' IJV disappears, 
as it ought to, from physical quantities such as the gravitational wave amplitude or the frequency and angular 
distributions of the emitted energy. 


3. RADIATION AMPLITUDE 

We proceed with the computation of the energy-momentum source of the gravitational radiation field. The 
gravitational wave source has two parts. One is the particle energy-momentum contribution, localized on the 
accelerated particle trajectories given in the previous section. The other is due to the non-linear self-interactions 
of the gravitational field spread over space-time. One should keep in mind that we are eventually interested in 
the computation of the emitted energy, given by dm It involves projection of the energy-momentum source 
on the polarization tensors and imposing the mass shell condition on the emitted radiation wave-vector. Thus, 
whenever convenient, we shall simplify the expressions for the Fourier transform of the energy-momentum source 
by imposing the on-shell condition k 2 = 0, as well as by projecting it on the two polarizations. 


A. Local source 


We start with the direct particle contribution to the source of radiation. We call it “local”, because, as 
mentioned above, it is localized on the particle trajectories. The first order term in the expansion of (12.61) is 


^Tnv{x) = e J do 2 l z^u v ) + 2 xu x h' x ^u v) - u^u^z • d) <5 4 (x -°z(cr)) , 

where z ^ is evaluated at a and h' is evaluated at °z IJ ’(o). Its Fourier transform is 

% v {k) = e ifc2(0) e [ doe dku)a fou (ll % + 2 xu x h! x{[1 u v) + i(k- x z) u M u„ 


Similarly for the primed particle with u replaced by u'. 
Introducing the momentum integrals 


1 = 


1 


S(qu') 5(ku — qu) e l ^ qb ' > 4 


(2ny 


d q, 


4 “ (2tt) 2 


6(qu r ) S(ku — qu) e l ( qb ' ) 


du di q > 


the first-order correction to the source becomes 2 

1 


% v = 2e 2 x 2 e lkz ^- n 

(. ku) 


( ku'I - + 2 U( M I„) 


\rpl _ lrri I 


(3.1) 


(3.2) 


(3.3) 


2 Terms, coming from the integration constants in (12.161) and (I2.17D . contain S(qu) and lead to the extra terms proportional to 
S(ku ) and S'(ku). With the on-shell condition k 2 = (ku)(ku') — k ^ = 0 the latter is equivalent to k^ = 0 and these terms do not 
contribute in the subsequent d?k— integration. 
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Note that the integrals I and Ip contain one massless Green’s function. This is in accordance with the fact 
that 4 Tp V , expressed through them, is the source of radiation from the colliding particles. I and Ip are computed 
in Appendix [Bl They are 




and upon substitution into (13.31) lead to 

= -2e 2 ^S kb ^ 2 

with apt! = (kb) UpU^ — 2 (ku) U(pb„y Similarly 


, (ku')m , .*'(&)*„. 
m "</■“"> + ~WuT + '-biku? 


U ) 1 

His 


^ = -2e 2 x 2 e~ i{kb)l2 


$(6) u\ 




(A:u)$(6) , , . ( E >, (&)cr^“ ^ 


2 (fax') 


■ u,,u ,, — l- 


b ( ku') 2 


(3.4) 


(3.5) 


(3.6) 


for the contribution of the primed particle, obtained from l T IJU by the substitution ► — b M ,u p -o- u ,p . 

Eventually, l Tp„ and 1 T / ( L/ will be contracted with the polarization vectors e± and e 2 , we will construct in the 
next section. They have zero time component and, therefore, satisfy e\ ■ u' = — e± ■ u and eg • 'u! = —e 2 • u. Thus, 
one may effectively replace in the energy momentum tensor u], by —Up when they are not contracted, to obtain 


and 


T — 

J /ll/ — 


-2 e a *V< fc W 2 


-®(b) UpU v + 


(ku')4>(b) .4>'(b)a$' 

2 (ku) U » Uv+l b(ku) 2 _ 


lrjif 

1 His 


—2e 2 x 2 e~ i ^ kb ^ 2 


-$( 6 ) UpU v 


(ku)$(b) .$'(&) affl' 

2 (ku') UpUl ' 1 b(ku') 2 


where a 


(«) = 
His — 


(kb) UpU„ + 2(ku') u^jb v y 


(3.7) 


(3.8) 


B. Non-local stress source 


The contribution to the source at second-order coming from the expansion of the Einstein tensor reads [I] 
Sfj,v(h) = h^’ p (h UP: x k^x p) h p (hp\^p T h v x,pp h\p p U hp^xp) 

~ h Xp ,phxp, v - ^ hp V d 2 h + i r]p V ^2 h xp d 2 hx p - hx p ,ah Xa ’ p + | h X p^h Xp ’ a j . 

It contains products of two first-order fields. Thus, it is not localized, hence its name “non-local”. It is also 
called “stress”, being part of the stress tensor of the gravitational field. 

Upon substitution of hp U and l z ,J ' (a) of the previous section in the above expression we obtain for the Fourier 
transform of Sp V 


Spv(k) = x 2 e 2 e l< ' kb ' )/2 (ku') 2 UpU„J + (ku) 2 u'pulJ + 4Jp„ + 4(ku') U(pJ v ) - 4(ku) u'^J v s ) + 
+2u(pu'^ (2 (kJ) — ( ku)(ku!) J — 2 Sp J^j 


in terms of the integrals 


Ja 


, 0 ) = 


(2n) 2 


S(qu') S(ku — qu) e 
q 2 (k-q) 2 


Qm -Qui 


d 4 q 


(l = 0,1, 2). We use the definition Sp J = ri^Jpv, while we have omitted the terms proportional to r]p V as well 
as the longitudinal ones proportional to kp or k u in anticipation of the fact that they will eventually vanish, 
when contracted with the radiation polarization tensors. Finally, as in the case of 4 Tp V one can effectively 
substitute u'p —> —Up to obtain: 


Sp„(k) = x 2 e 2 e ^ kb)/2 


^4Sp J — 4 (kJ) + [(ku') + (ku)] 2 j) UpU v + 4Jp V + 4 [(ku') + (ku)] U(pJ„) . (3.9) 
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Note that contain the product of two graviton Green’s functions, which signals the fact that S^ is due 
to radiation from “internal graviton lines” in a Feynman graph language, through the cubic graviton interaction 
terms. It will be explicitly demonstrated below that in the zero frequency limit the contribution of in the 
emitted radiation is negligible, as argued in ||. Nevertheless, it will become clear that it contributes significantly 
at high frequencies and, as will be shown next, it plays an important role in the cancellation of the r o dependence 
in physical quantities. 


C. Cancellation of the arbitrary scale r q 


As anticipated, in this subsection we will demonstrate explicitly that the arbitrary scale ro disappears from the 
final expression of the total contribution to the source l T^ liV + l X)(„ + S^ v of the gravitational radiation. As will 
become clear below, the local and stress parts of the source each depends on r o, but their sum is ro—independent 
and finite. According to their expressions in (13.71) and (13.81) . l T IJL , and depend on r 0 through 4>(&), while 
Sfn, depends on ro through terms proportional to AT_ i(£) (with no extra factors £) in the expressions of J, 
and J ^ v , evaluated in Appendix [B] All these unphysical terms will be shown to cancel out and will end up with 
expressions (13.241) and (13.251) for the total energy-momentum source for the two polarizations separately 3 . 

We proceed in steps: 

1. Split Sp V = S 1 ^ + S^ v with 4 


S£„ = 4^ 2 e 2 e i(feb) / 2 (Sp J - [k J)) 


U^U v 


S%, = x 2 e 2 e ? m/2 


1 [{ku') + ( ku)] 2 J u,j,u u + 4 J M1/ -f 4 [( ku ') + (ku)] u^J l 
Using (IeTgI) and (IB.8j) . becomes 


V) 


Q 1 


-x 2 e 2 §(b) e j(fcb)/2 (e _i(fcb) + l) = -2 x 2 e 2 $(b) cos 
2. Similarly, it is convenient to split the local source l T llv + (13.7113.81) 

e i(feb)/2 + e -i(kb )/2 


y . 


T t l = e 2 x 2 


T n = -e 2 

^ fJIS ° 


rnlll _ 


vl 

2 


as: 
k • b 

4>(6) u^Uv = 2 x 2 e 2 $(6) cos —— 


i(kb)/2 ( ku ') —i(kb)/2 ( ku ) 

(ku) (ku 1 ). 




,&(b) 


(3.10) 


(3.11) 


(3.12) 


Thus, + S 1 ^ = 0. 

3. The remaining stress contribution S 1 ^ is a linear combination of J, J^ and J^ , which have been computed 
in Appendix [Bj Taking, as above, into account the fact that they will eventually be contracted with the 
polarization vectors and that we shall set k 2 = 0 in the integral for the radiation energy and momentum we are 
interested in, they are 5 : 


1 

J= — / dxe-W^K - x (k±by/x(l - x)) , 


T eS — _ / p -i{kb)x 

^ ~8n 


NfK^iO+i^KoiC) 


Nf = — . 
^ 2 L 


x(ku') + (1 — x)(ku ) 


Jf, = y / dx 


^NfNfK.k C) + (2 iNtfb v) - u^KoiO - ^yA'i(C) 


3 The reader, who is not interested in the details, may go directly to these formulae for the total source. 

4 The integrals k • J and Sp J are singled-out, because they can be computed exactly. See Appendix [B] 

5 Note that we use non-standard symbols for the modified Bessel functions, namely K u (z) = K u (z) z v . In this notation the 
differentiation rule reads K^z) = —zK u —i for any is, while the zero-argument limit is AT(0) = 2 U ~ 1 T(i/) for is > 0. 





























Having anticipated that the dange rous term s for divergence and r $—dependence are the ones which contain the 
integral of K_ i(f) with ( = k±b^/x(l — x), since according to Appendix [Bl lead to $( 6) 6 it is natural to treat 
separately the terms in S Ik which contain K-\, from the ones which contain Kq or K\. Thus, in a suggestive 
notation, we split: S 11 = S^) 1 ' + with 


S'- 1 ) = 2 Gb 2 e 2 e* kb) / 2 / \ [( ku ') + ( ku )] 2 u M u„ + 4 NfNf + 4 [(ku') + ( ku )] u ( „N°f 


S,'?/ 1 ) = SGe 2 e i(kb)/2 dxe~ i{kb)x 


- UfjUv + i [( ku ') + (ku)] u^b^K 0 (0 - ^^A'i(C) 


■ (3.13) 


Substituting the explicit form of iV® ff and simplifying, we obtain 


S^ 1 ) = 2 G 6 2 e 2 e i(fcb) / 2 / dxe~^ kb)x K-^C,) [(1 - x)(ku') + x(ku) 


U^Uv . 


(3.14) 


4. Consider, next X/JJ ,. Using the formulae derived in Appendix [Bl i.e. 


1 

e~^ kb ^(b) = ±-Jdx e" i(fcb)x [x(2x - 1 ) b 2 k i^-r(C) + 2 (l - ix(kb)^K 0 ( C) 


1 

$(&) = -L J dx e- i{kb)x (x — l)(2x — 1 ) b 2 k 2 ±K-\(Q + 2 ^1 — i(x — l)(fcfe)^Ao(C) , (3.15) 


T,{{. takes the form 


T" = - 2 Gb 2 e 2 e l{kb)l2 / dx e ~ i(kb)x 


(x — 1 )(ku') 2 + x(ku) 2 (2x — l)A"_i(^)+ 


(3.16) 


+ 2 


(ku') 2 ^1 — i(x — 1 )(kb)j + (ku) 2 ^1 — ix(kb)^ 


Ko(0 


k^b 2 


u^u v . 


5. Thus, the sum 


S*" 1 ) + T 11 = - 2 Gb 2 e 2 e i(kb)/2 / dx e" i(feb)a 


(ku') 2 (l — i(x — l)(fc&)^ + ( ku ) 2 ^ 


1 — ix 


Ko(() 

J k\b 2 


(ku') + (ku) x(l — x) K_i(() 




(3.17) 


and, using the identity z 2 K-i(z) = Ki(z), we obtain the explicitly finite expression 

l 


S'!” 1 ) -L J'U — i(fcb)/2 

a fiv > ± fiv — jZ2 e Ufj,u v 


dxe~ i{kb)x 


1 2 


(ku') + (ku) Ki(()— 


(3.18) 


-2 


(ku') 2 ^1 — i(x — 1 )(kb)^j + (ku) 2 ^1 — ix(kb)j Ko(C) 


with no A'_i(<C). All divergent and ro~dependent tei'ms have cancelled. 


The integral containing the hatted Macdonald of index —1, which near x = 0,1 behaves as K—\ \k±by/x{l — x)j ~ [x(l — *)]-\ 
diverges logarithmically at both ends of the integration region. In Appendix |B] it is shown that this logarithmic behavior is 
related to the one of <f> (Eqns. JB.6I [B~8l) ^ . Alternatively, one could regularize these divergent integrals by shifting the index of 
all Macdonald functions by 0 < e 1, which makes all integrations convergent, and take the limit e —> 0 in the very end of 
the computation. 
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D. The total amplitude 


Therefore the total effective radiation amplitude reads 

T — 9(0,1) 1 o(-l) 4. ' T 11 1 'T 

1 /J.V — 1 -L fjn, 1 -L fju, 

• From its definition in (13.131) . Sfm’ 1 '* takes the form 


hi 


(3.19) 


Sf& 1) = 8Ge 2 e l(fcfc)/2 / dxe~ i{kb)x 


-u^u v K 0 (Q - ^^Ad(C) + i (1 - x)(ku') + x(ku) u^b^K 0 (Q 


• Using formulae 

p -i(kb) ®'( b ) 


4nb 2 


1 

J dx e-^Uz x(2x - 1) b 2 k 2 ± K 0 (() + 2 (l - ix(kb)) A'i(C) 


$'( 6 ) 


47 xb 2 


1 

Jdxe~^ x (x — l)(2x — 1) b 2 k^_Ko(C) + 2 ^1 — i(x — l)(kb)) Ki(() , (3.20) 


l T^l in (13.121) is also written as a sum of two integrals over x, one containing K 0 (Q and the other Ki{(). 

• Collecting terms with integrand proportional to K 0 and Ad we write the total energy-momentum source 
t™ in the form 


r^{k) = 2 Ge 2 e l(fcb) / 2 / dxe~ i{kb)x \ Y°(k, x) AT 0 (C) + YUk, x) Ab(C) 


(3.21) 


where 


F u = - 4 


ki L 


( ku ') 2 + (kuY 


4 i 


. (kb)x( 1 — x) 


k\ 


+ 8fx(l — x) ( ku') + (ku) 




Y ^=(ki 


(ku') + (Aru) 


4z 


, m 


k\b 2 


(ku') 2 — (kuY 


(ku') 2 — (ku) 2 

(kb) 2 r 


U^Uv- 


-4 


L 


(1 — x)(ku') 2 + x(ku)' 




-4 


b^bu 8 


b 2 k 2 ± b 2 


(i (ku') + (ku) + (kb) (x — l)(ku') + x(ku) ) u^bu) ■ 


(3.22) 

(3.23) 


• Defining in the center-of-mass frame the radiation wave-vector by 

= w(l,n) = uj ( 1 , sin d cos gs, sin d sin gs, cos d) 

and contracting (13.211) with the two polarizations (see Appendix 0, we obtain the final (finite) expressions for 
the source of the gravitational radiation separately for the two polarizations 7 

1 


r + (k) = r M! y(fc)e^ = 


16 Ge" 

~7T 


dxe~ i{kb)x 


-K 0 (C) +sin 2 ^ATi(C) 


and 


r x (k) = r A1! y(fc)e x l/ = — sings [ dxe d kb ) x 


V2 


X * 2 ® . Al(C) + (2* - l)cosgsAd(C) 

ub sin v 


(3.24) 


(3.25) 


where ( = uib sindy 7 x(l — x). 

To summarize: The only approximation made so far is the restriction to the first order corrections of the 
gravitational field. The leading non-linear terms were taken into account. To this order, the total source r 
H3.2H) and the separate sources of the two polarizations \3.24\) and Iff. 25\) have been expressed as finite integrals 
over a parameter x € [0,1]. 


7 Up to an overall phase e l (k&)/ 2 in both r_|_ and r x , since it does not contribute to the energy. 
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Figure 1: The characteristic angular and frequency regimes. 


4. CHARACTERISTICS OF THE EMITTED RADIATION 


We turn next to the computation of the emitted radiation frequency spectrum and of the total emitted energy. 
They are obtained from 


d Traci 
dco dfl 



J2\ t -p I 2 * 


(4.1) 


summed over the two polarizations. 

It will be convenient in the sequel to treat separately the six angular and frequency regimes shown in Fig. [I] 


A. Zero-frequency limit — Regimes I and II 


In the low-frequency regime (co —» 0) the amplitude r x dominates and has the form 


_ 16V2iGE 2 siny f , i(kb)x 
x cobsind I d 


K 2 (C) - K x { C) 


(4.2) 


while t + is finite and gives subleading contribution to (14.11) . Note that in this limit (/\i(0) = 1, K 2 {0) = 2, 
e -i(kb)x _ the a;—integration is trivial and gives dE lad /dco = (2 8 G 3 E 4 /nb 2 ) J did/ sind, which diverges and 
implies that our formulae are not valid for id close to zero. 

We cannot trust our formulae in that regime and should repair them. A quick way to do it, is to impose a 
small-angle cut-off id = id cr on the id— integration, so as to obtain for d.E ra( j/da;| w= o the value computed quantum 
mechanically in [s|, Q, namely 


dE w \ 

^ ) u ,=0 




Ga 2 s 

ti 


4 Eb 

In — =- a 3 ln(2 /a), 

or tt 


which indeed agrees with our expression 


s = 4E 2 , 


tt/2 


( dE lad ' 

\ 2 9 G 3 E 4 

r did 

V dw , 

L=o nb2 J 

1 sind 


l?cr 


Eb 3 idcr 

— a° lnctg — 
7r Z 


(4.3) 


(4.4) 


for 


id cr = a « 1 . 


( 4 . 5 ) 
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Thus, our result for the low frequency radiation emitted in a collision with large s, fixed t = —(s/2) (1—cosa) ~ 
—sa 2 /4 and G\t\ ^ 0(1), agrees with the quantum computation of Weinberg, apart from a tiny emission angle 
id ^ 0(a) -Cl in the forward direction 8 . Furthermore, as will be shown next, at low frequencies oj and for 
d ^ a, the leading contribution to our classical amplitude, dominated in this regime by Tpf \3.1 6 2\) . is identical 
to the one obtained in J 6], after it is generalized (see below) to b ^ 0 9 . 

Indeed, following the notation of we write for the energy-momentum source of the 2 —» 2 scattering process 
we are studying for arbitrary, a priori, scattering angle a and impact parameter b: 


2 9 ~ ~ 

pft pis pfi pis 

T‘ “"(x, t) = J2 -IT 1 - <* 3 ( x - v » 4 T b/ 2 ) 0 H) + Y ^( x 

_i En _i E n 


A pi-L piS 

= 'E^TT LS3 (*- V » tz Fb/2) + E 


( P^p 

I n ? 


n r3 


n =1 


V En 


6 6 (x-v n t- Fb/2) 


T b/2) 6(t) 

PM pv \ 

-^i 3 (x-v„iTb/ 2 ) )d(t) 
£sn J 


where = 1,2, are the initial (final) particle momenta, with P n = P n (P n ,a). Its Fourier transform is 

2 2 


f^(k,u>) = Y 2 tt e ±l(kb)/2 P£P" S(k -P n ) + Y e ±i(/cb)/2 


f pupv 2 
V E n oj — kv„ 


pn p v » 

E n w - kv„ 


where = (w, k) is the radiation wave-vector. The terms n = 1 (n = 2) in the sums, are multiplied by e+d kb )/ 2 

( e —i(fcb)/2), respectively. 

The first sum, proportional to delta-functions, corresponds to no scattering and does not contribute to 
radiation. Thus, we end-up effectively with 


J / pn p i/ 

T^(k) = i Y e ±i(fe&)/2 I —— 


n—1 


V E n 


- kv„ 


pn pv | 

n n _ 

u - kv„ 


(4.6) 


To leading order in our approximation the scattering process, we are dealing with, is elastic with E n = 
E n = E. Furthermore, write for the incoming particles = Eu% = £1(1,0, 0,±1) and for the outgoing ones 

Pf) = Euif = E(u% + Hff) = E(u^ =F a 6 M ) = P)) + x Pn, substitute into , and expand in powers of a using 
the fact that for d > a one has |k ■ 3 P%\ <C \k -Pf\, to obtain 


T%(k) = iY e 


zi{kb)/2 ( 2Pn' J ' 1 Pn ) _ M^n) ^ v \ 0(0?) 

\ (kPn) (kP n ) 2 n U ' + ( ’ 


(4.7) 


or, finally, making use of the above definitions, 




^±i(kb)/2 


(kb )« - 2 (ku n ) u^b u) _ 8iGE 2 
(ku n ) 2 b 2 


(“) (“') 

J(kb)/2 a ni' _ -i(kb )/2 ' 

(ku) 2 (ku') 2 


(4.8) 


This is identical to Tj£ in (13.121) . after we bring the latter to its original form (13.61) by the substitution ajfi) —»■ 
atf\ Q.E.D. 

Summary: In regime II we use our formula, which is identical to Weinberg’s. We extend the use of Wein¬ 
berg’s formula in regime I as well. The total emitted energy in I is known to be of 0(a 3 E). Given that 
(dE ra< i/duj) u -o ~ (dE ISL d/dui) u= i/i ) , the contribution of II to the radiation efficiency is estimated by multiplying 
( 1 4-4V by the frequency range 1/b and dividing by the initial energy 2 E. The result is 

eiji ~ a 3 ln(l/a). (4-9) 


8 It should be pointed out that our classical computation reproduces the quantum results of Weinberg for soft graviton emission 

for G\t\ < 0(1), i.e. for (E/Mpi)(rs/b) 1, or in Weinberg’s notation for B 1. 

9 It is not surprising that the quantum and the classical results agree for the emitted energy of low frequency. As we argued, in 
the low-frequency regime, dE ra d/dw is dominated by the local source. The contribution of stress in this regime is negligible. As 
a result, the radiated gravitons are expected to be produced in a coherent state m. and the corresponding expectation value of 
the quantum field to satisfy the classical field equations. 
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B. Regime VI 


Here we consider the Regime VI (w > l/(a&) = 1/2 rs and $ > a). Contrary to regimes I and II, in regime 
VI as well as in IV, which will be discussed in the next subsection, the contributions to radiation of the local 
term and the stress are equally important. 

The cross-amplitude (13.251) after an integration by parts of the second term reads: 


32 iGe 2 

t x (k) = - -y=— sin ip 


' dx e -i(kb)x ( K 2 {Q^ 2 <P-Ki{Q \ _ 4 COS if e _ i(fcb)/2 sin kb 


with a = ujb sin d and together the two amplitudes are written as 

16Ge 2 - 


T + (k) = 


V2 


,,, 32 iGe 2 . 

Tx{k) =- -j=— sinyj 


— Lq + L\ sin 2 y> 

L 2 sin 2 y> — Li 4 cos ip p -i( kb )/ 2 . q kb 


where 


l 

L m (a,ip) = f e mx cos v k m (a\/ x(l - x)^j dx 


is defined and studied in Appendix [C] 

Substituting the large-a expansion to leading order 


L r n — 


2 m+2 T(m+ 1) _ iacoBv/2 co ~ acosip 


(4.10) 


we obtain 


Thus 


64Ge 2 acosip 64 iGe 2 . . acosip 

r + (k) ps- -= — cos 2 <pcos —-— , r x (k) ps- sm2 ip sin —- — 


y/2 a 2 


dE lad 2 10 G 3 A 4 u 2 


V^a 2 


2 r,/a cos ip\ . o . r,/acos(p\ 
cos 2 2 y;cos 2 ^———J + sin 2 2 y>sm“ ——— J 


Integrate over ip using the formulae 1 A. 71) to obtain for a > 0( 1) 


dE IBd _ 2 10 G 3 E 4 uj 2 
duidd it a 4 


(l-™) J 2 (a) 
a V cl 


. 2(ab) 3 Eu; 2 . 

sin v ~- j sin v , 


(4.11) 


(4.12) 


(4.13) 


from which one can obtain an estimate for the frequency distribution of the emitted radiation in regime VI by 
integrating over d £ (a, n — a), namely 


dEZ aE 1 


ui > 1/ab , 


(4.14) 


as well as an estimate for the emitted energy and the corresponding efficiency in regime VI, by integrating also 
over ui £ (1 /ab, oo), 


E, 


VI 

rad 


l E and 


€yi ~ a . 


(4.15) 


C. Regimes II+IV+VI 

Finally, we shall discuss the characteristics of the radiation corresponding to the union of regimes II+IV+VI. 
For that, it is convenient to transform the amplitudes (13.241) and (13.251) to the form 


T +( a ; if) = 


16 Gee' 


dx e iax cos f 


-K 0 {£) +sin 2 ^Ad(C) 


(4.16) 


Tx (a, <fi) = ~ 


16 i Gee' 


sin p / dxe 1 ' 


2 a;(1 — x)aK 0 (Q + ( —i{2x — 1) cosy; + - ) A'i(C) 


o 
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Representing Macdonalds as 


A- o( o=/ exp( ~i^d I) ^, 


0 


VW+T 


Ki(0 = C 2 


exp (-C Vt 2 + 1 ) 2 


\ft T n 


t dt , 


(4.17) 


we square and sum up the two polarizations: 

l l 


OO CO 


/ p p p — i) x^ j 4 

dx J dx ' J dt J dt ' Vt 2 + Wt' 2 + 1 'Yh Mk ( x V, x ', t ',v) 

n n n n k—0 


0 0 


0 0 


A(a;, t, x', t', p) = \/x(l — x)\/t 2 + 1 + \/x'(l — x')\Jt' 2 + 1 — i(x — a/) cos 




(4.18) 


with 


M 0 = l 

Mi = 0 

Af 2 = a;(l — s) sin 2 y> —t 2 + 2 a/(l — a/)(l + 2 t 2 + t 2 t' 2 ) + {a 


x', t 


M 3 = 2ixx'(l — x)(l — x') t 12 sin 2 <^cosv? 


2x' — 1 + 2x't 2 


— {a: •<— > x ', f 


*'} 


M 4 = xa/(l — a;)(l — x') t 2 t' 2 sin 2 </? sin 2 </? + ( 2x — l)(2x' — 1 ) cos 2 p 


(4.19) 


One may substitute into (sm and integrate numerically over all variables apart from ui to obtain dE rad /dui. 
The result is shown in Fig. [2] On the other hand, one may change variables and integrate over a instead of to 
to obtain: 


dE, 


rad 


dd 


2 6 G 3 E i 1 
n 2 b 3 sin 2 $ 


4 i i oo 

^^(fc + 2 )! f dx f dx' f 
*=o o o n 


dt 


dt' 


V^Ti J Vt' 2 + 1 


dp 


M k 

\k +3 


The convergence of this 5-variable integral follows from the behavior of the integrand at large-a, while its 
numerical value is expected to be of 0(1), since the integrand does not contain any small or large parameter. 
Indeed, numerical integration (including the factorial in front) leads to 


oo oo 2n 


r l = ^2(k + ‘2.)\ J dx J dx' J dt J dt' J dp-^-^= 

k—0 0 0 n n n 


M k 


o o 


2 + lv't' 2 + 1 


89.9. 


Thus the angular distribution reads 


Integration over £ (a, n — a) gives 


dE rad _ ija 3 E 1 
dti 


E ra d — 


87T 2 sin 2 $ 
r] a 2 E 


47T 2 


and for the efficiency 


■ = w = Lli « ! 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


Finally, the ip —distribution of the emitted radiation is shown in Fig.[3l according to which most energy is 
emitted perpendicular to the scattering plane. 

It is instructive to study regime IV in a little more detail. For that, let us split regime IV into IVa and IVb, 
as shown in Fig. [I] according to a > 1 and a < 1, respectively. Inside the regime IVa the amplitude is damped 
as in regime VI. However, near the left border of regime IVb (with 1/b < to <C 1/ab) one may expand the 
amplitudes in powers of a and obtain, as in regime II: 


r~r y ~ - 


16-\/2 iGE 2 sin(/j 


ujb sintl 


(4.24) 
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Figure 2: The frequency distribution in the combined regimes II+IV+VI for a = 0.01 with two different choices of axes’ 
labelings. The regime II in the Figure on the left is compressed in the interval (0, a). The slope in regime VI of the right 
plot corresponds to dE la , d /doJ ~ 1/cu 2 . 



Figure 3: The ip— distribution for G = b = 1 and a = 0.01. 


Upon integration over regime IVb, i.e. for a < i9 < i9 max = arcsin(l/w6) one obtains 


f ^ 1-Vad \ 


Eb tg(i? max /2) a 3 Eb 
in 


In- 


2a- 


V du J i/6<u,«i/a6 n tg(a/2) tt U>b + Vu 2 b 2 - 1 

Thus, for 1/& < w <C l/a6 one may approximate dE lad /du} by 

/ dE lad \ ^ a 3 Eb ^ 2 

V dto A/6<u;«l/afc _ * Otub ' 

On the other hand, from the known behavior near 1/ab from regime VI, we know that 

dE r a d 


du) 


lEb. 


I cub 

So, a natural interpolation of dE rad /dui between the values 1/6 and l/a6 is 

(dE^ aA \ _ a 3 Eb . 2 

V du} )i/b<u<i/*b n ^ n aujb 


(4.25) 
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It may be shown numerically that £(w 6 ) is a slowly varying function of 0(1) in the regime 1/6 < ui < 1/ab, so 
that one can simply write instead 


/ dE r a.d \ 

V diO ) 1/6< W <1/ Q i, 


a 3 Eb In 


auib 


(4.26) 


It should be pointed out here that the integral of dE la d/dco over ui receives most of its contribution from 
frequencies in the neighborhood of l/a 6 in both regimes IV and VI. Thus, one can say that the characteristic 
frequency of the emitted radiation is around 0(l/rs). 


D. Comparison to previous work 


Let us briefly compare the results of the present paper with our previous calculation of massive particle 
collisions and to the previous literature as well. In relation to our work, we should point out that it is not 
straightforward to compare with [lj, simply because they are were based on different assumptions and approxi¬ 
mations. In particular, in the massive case we worked in the lab-frame and assumed weak-field approximation, 
restricting 7 cm to 7 cm < 1 /a, which is not consistent for fixed a with the 7 cm —> oo limit, necessary to connect 
with the massless case. On the other hand, in the massless case the computation is characterized by 7 cm = oo, 
is performed in the center-of-mass frame, used the scattering angle a as an angular integration cut-off, and used 
Weinberg’s results in regime I. So, before it is appropriate to compare with the massive case, we have to redo 
the latter following the same methodology and approximations. This, as well as comparison with other papers 
in the literature is the content of this subsection. 

In [l] we computed the gravitational radiation efficiency in ultra relativistic massive particle collisions in D 
dimensions and in the weak field limit, i.e. for 1 <C 7 cm •C 1/a. Specifically, in D = 4 we obtained e ~ or 7 C m, 
with the emitted radiation having characteristic frequency w cm ~ 7cm/6 and emission angle d ~ l/7 cm , where 
6 is the impact parameter. Thus, the efficiency is finite in the above kinematical regime and agrees with the 
results of [lj| and jl6j |. 

However, to compare with the massless case we have to improve the massive case computation in two ways: 

(i) discard the weak-field condition as it is done here and more generally in the literature, e.g. ®mm , and 

(ii) use the scattering angle a as a characteristic angle to cut-off angular integrations. As a consequence of (i) 

the condition 1 <C 7 cm •C 1/a is replaced by 7 cm 1 and 1/a 1 , independently; so that for 7 cm —> oo one 

would naively conclude that e ~ a 3 7 cm —> oo, in disagreement with the present massless result e ~ a 2 . But, as 
we will argue next, the formula e ~ a 3 7 cm is not valid in this case and the present result from regimes I-IV and 
VI, defined by a , 1/6 and l/a 6 as in Fig. [Q is actually in agreement with the limit 7 cm —> oo of the contribution 
from those regimes in the massive case. 

Let us start with regime IV a + VI. Due to destructive interference [l], the total radiation amplitude is 


T(w,'i?) 


Gm 2 7/m 
(ujbsind) 2 ’ 


(4.27) 


which upon integration over this region gives e ~ a 2 + O (l/ 7 cm) • 

The contribution from regimes II + IVb can be computed as in the massless case. That is, starting with the 
observation that in that regime the direct amplitude T dominates over the stress S , we take for T the expression 
0 , Eqn. (3.17)] 10 




iGm 2 7 ( 


6 2 


2 

■£ 

n =1 


iviy) 


±e : 


±i(kb) /2 (fc 6 ) « - 2{ku n )ujb 
( ku n ) 2 


AT 


( ku n )b 
SOori 


(4.28) 


and integrate over the above regime to obtain again e ~ a 2 , plus subleading corrections of order l/ 7 cm ■ 

Next, one can, following the approach discussed in subsection 14 Al above, use instead of (14.281) the expression 
XT/ given in (14. 61) . valid for both massive and massless colliding particles. This allows us to integrate over all 
angles and upon integration over I+III to obtain a gain to e ~ a 2 . Furthermore, this approach eliminates the 
disagreement of 0, Eqn. (4.15)] as well as of Smarr |14l Eqn. (2.13)] in the expression for dE la d/du> in the zero 


10 In this expression the true 4-velocity U^ normalized by U 2 = 1, is defined as U^ = 7 cm w M with 2 = (1,0, 0, ±v C m) and 
^cm ~ 1 — 7cm • 
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frequency limit with the corresponding result of Weinberg. Specifically, the non-sensical for 7 cm —> oo formula 
(dE lad /du!) u=0 — (Gsa 2 /ir) ln( 47 2 m ) gets modified to (dE^/cko)^^ — (Gsa 2 /n) ln(s/|t|) given in (POl) . 

Finally, in regime V the direct amplitude and the stress are equally important, both for massive as well as 
massless colliding particles. This was important in deriving destructive interference, but on the other hand it 
does not allow us to replace the source amplitude with Weinberg’s formula for T. To proceed along these lines, 
one has to “guess” the correct corresponding modification of S, consistent with the conservation requirements 
for the total energy momentum tensor. This problem is currently under investigation. 

Thus, indeed, the results of the present paper can be obtained as the massless limit of the corresponding 
conclusions of the massive case in their common region of validity. 


5. CONCLUSIONS - DISCUSSION 

Using the same approach as in [l|, based on standard GR, with the leading non-linear gravity effects taken 
into account, we studied collisions of massless particles and computed the gravitational energy of arbitrary 
frequency, which is emitted outside the cone of angle a = 2 r$/b <C 1 in the forward and backward directions. 
The value e ~ 1.14 a 2 was obtained for the radiation efficiency, with characteristic frequency w ~ 1/rg. In fact, 
this value represents a lower bound of the efficiency, since it does not include the energy emitted inside that 
cone. The frequency distribution of radiation in the characteristic angle-frequency regimes is shown in Fig.[T] 
Our method allowed to study the zero frequency regime and showed that the values dE lad /du> at u> = 0 and 1/6 
are of the same order. 

We would like to point out that our results about regime IV agree with Gruzinov and Veneziano 0 - Fur¬ 
thermore, our work provides information about the very low frequency regime II, in which, strictly speaking, the 
method of [l0| cannot be applied. Finally, at frequencies ui < l/a 6 we obtain e ~ a 2 , i.e. the same dependence 
as in m- A possible difference in the coefficients is expected, since the non-linear stress contribution around 
ui ~ l/a 6 is comparable to the contribution of the direct linearized emission from the colliding particles. How¬ 
ever, as mentioned several times, our approach cannot yet deal reliably with regime V, which according to [10] 
gives dominant contribution in e by an extra factor In 1/a. In particular, we cannot yet confirm the presence of 
any other characteristic frequency, such as e.g. l/a 3 6 , or characteristic emission angle smaller than a jiol ]. We 
hope to return to these issues with a better understanding of regime V in the near future. 
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Appendix A: Conventions 


Our convention for the Minkowski metric is r = diag(+l, —1, —1, —1). 

Fourier transforms are defined as: 

/(*) = 7^4 / /( fc ) e~ i(kx) d 4 k , f(k) = J f(x) e^d 4 x = E[f(x)](k) . (A.l) 

We use the following symmetrization notation for tensorial indices: a^b^ = [a M 6 „ + a u b^]/2. 


Polarization tensors. Given the radiation null wave-vector k p and the null velocities u p and u' p , define 
the polarization vectors ei and 62 by: 


k _l 


(fcu)u' M — {ku')u p ku — ku' 


ku + ku' 


k " 


2 kj_ 


e puXp u v u' x k 


P > 


(A.2) 


where k\ = (, ku)(ku ') and e 0123 = 1. They satisfy 

ei • e 2 = k ■ e\ = k ■ &2 = 0 and e 2 • u = e 2 • u = e 2 • k = 0 . (A.3) 

In the center of mass frame of the collision under study we have chosen u p = (1, 0,0,1), u' p = (1,0,0, —1) and 
6 ^ = (0,6,0,0). In addition, we have defined k p = w(l,n) = w(l,sin $ cost/?,sin t? sin:/?,cost?). Using these, we 
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obtain 


e(* = ( 0 , dn/d'd ), eit = -- ( 0 , dn/dip) 

sm v 


and 


ei ■ u = sin i9 , e\- u! = — sin i9 , e\ ■ b = — b cos $ cos <p , 

Finally, using ei and e 2 , we define the polarization tensors as 

£x" = V2e^e'J , r)^ = 0, 


LI V LI V 

fj,v _ e l e l — e 2 e 2 

+ " V2 


e 2 ■ b = — 6 sint/ 3 . 


V = +,x. 


(A.4) 


(A.5) 


(A. 6 ) 


Useful integrals. The following integrals are used in the main text: 


2 tt 


1 Ji(a) 

= 71- | -- 

2 a 


/ . 2 . 2 facos(p\ 

/ sm ipsm ^—-—J dtp 

o 

2 tt 

f . 4 • 2 f aCOS( P\ i of 1 M a ) 

j emv sm (-j-)<f* = 3,r (j- — 


j ' J ' j = 7T Q - J 0 {d) + 


0 

2 tt 


4 . 2 / a COS if \ 

cos (/ 9 sm ^—-—J dtp 


1 “ 72 J 2 ( 


(A.7) 


Appendix B: Computation of momentum integrals 

It is convenient for the computation of the integrals below to define the four vector by: 


b^ = 


1 

2 L 


(6m) m ,m + (6m') m m 


+ 6 /. 


Local integrals. The basic scalar integral is 
1 


/ = 


(2tt)2 


<S(gu') <5(/cw - gu) e~^ ^ _ <^i{ b u')(ku)/2 f p i q b 


Q 


2(2tt) 2 


( B-1 ) 


with b = \J—b ■ b and d>(r) = —(1/27r) ln(r/ro). 

The local vector momentum integral is defined as 


Jm ~ (2tt) 2 


6{qu')6(ku-qu)e ^ 4 , <91 jf 6 U 'lffcul /2 

q^ a q — 1 aun - e 


(fcn)$( 6 ) , .$'(&) , 

4 1 2b M 


g■“ ■ _ <96^ 

The corresponding integrals for the primed particle are obtained by the substitution m m o m , ai , b M 
Stress integrals. The covariant stress integrals are defined by 


J = 


1 f S(qu') 5(ku — gu) e 


(27T) 2 


g 2 (fc-g ) 2 


,4 J =i 9J T 

d q, Jn — 1 > dfj, v — 


d 2 J 


^ db^db" ' 


(B.2) 

- 6 ^. 

(B.3) 


• We start with the computation of J. We will then compute the other stress integrals as derivatives of J, 
following the same procedure as above, but to avoid unnecessary complications we will give the final expressions 
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only for the special choice 6^ = (0,1,0, 0) made in the main text. We obtain 11 

l 


J = 


1 


v'qfa 


7T OO 


87T 2 J q 2 (k — q) 2 


d 2 q = 


1 

2-7T 2 


dx e ~ i{kb)x J de J dq 


r^iqb cos 6 


4-7T 


dx e- i(kh ^ x J dq — 


qJ 0 (qb) 


o o 
i 


[q 2 + k 2 a:(l — x )] 2 


q = 


2 + k 2 x(l — a;)] 2 87 t 


J dxe-^^K- 1 (k±by/x(l - xfj 


where k± = |k| = \J ( ku)[ku') is the magnitude of the two-dimensional transverse part of 
Also 12 , 


^ = — / dx e -^ x 


NuK-xiQ+i-gKoiC) 


N 11 =xk^-- v 


x(ku') — (1 — x)(ku) , (B.4) 


where £ = k±b^/x(l — x). 

• Next consider (kJ). With the on-shell condition k 2 = 0 it can be represented as 

t - j - / g w d\=\j ^ 5 J£l - g ^ _ (t _ qf] 


1 

2 L 


e ~i(kb)j* _ j _ _ —i(fcb) _ ^ 

J 2 L 


$(&) 


1 _ e ~i( kb ) 


On the other hand, the integral form (IB. 41) of becomes 

l 

i k\b 2 K-i(jc±b\/x(l — x + «(/c6) A'o (k±by/:x{\ — x)^j 

whose integrand is the derivative of —e~ l ( kb ) x Kq (k±b\/ x{\ — x)^j . 
k • J = -L e- i(feb)x A:o(fcj.fo\/a;(l-x)) 


k-J = — I dx e ~^ kb)x 
SnJ 
0 


Hence, 


° = lim A> 2 (0) 
i d —^+o 87 r 


1 _ e -*( fcb ) _ J _ e ~i(kb) 


• Following similar steps, we find for the tensorial integral 
1 


= — / dx e-«‘ 


^AyWA-itC) + (2*JV (#i & I , ) - + u ( m <)) tf 0 (C) - ^ A'i(C) 


(B.5) 


(B.6) 


(B.7) 


Its trace Sp J = r) ,lu J,, v is given by 


Sp J = — dx e~ l(kb)x 

87T ./ 


-xWklKMO + 2ix(kb)K 0 (() + KM) - 2K 0 (0 


11 Use was made of the formulae 
r ±iz cos e oin 2 k PJQ _ 2 fc r(/c + l/2)y / 7T 

b 

12 A useful formula is J° =-J z = J. 

2 


TV 

/■ 


JfcW 1 


oc 

/ 


^ [ y 2 + 0 2]m - 2—lr(m) +l( b l a D- 
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Make use of the Macdonald-function identity K\(z) = z 2 K_i(z) to obtain 

1 

— x(2x — 1) b 2 k\K- 1 (£) + 2 (ix{kb) — l)/vo(</) 


Sp J = — dx e~ i(kb)x 

87T ./ 


o 


The integrand is the total derivative of —2*e z ^ kb ^ x K Q (jt±b^yx( 1 — a;)) . Hence 

Sp J = ^ e-^ kb)x K 0 (k±b^/x(l - x)^j ° = - e - i(fe6) ® = e~ i{ - kb) I. 


(B.8) 


Appendix C: Asymptotic expansion of higher-frequency integrals 

We are interested here in the expansion of the typical integrals 

l 

L'tri = J Z~ lCX Km (ay/x(l x)^j dx 
0 

(where c = k ■ b = — a cosip) in powers of 1/a for a^$> 1. The presence of the \/ x(l — x) in the argument of 
the K m , prevents the direct use of an asymptotic expansion of the Macdonald function itself. Instead, we shall 
proceed as follows: we represent 


— 


._ OO 1 

\fl (|) J dyK m+1/2 {| \/y 2 + l) J e~ lcx ch (ay (x - 1/2)) dx , 


where K m (z) = z m K m (z). Making use of 

l 


: cli^s (x — 1/2)) dx 


2 e~ ic / 2 


. s c , s . c 
s sh - cos - + c ch - sin — 
2 2 2 2 


2 \ 2 
S ^ C , 


one obtains 


Lm. — 


2e —ic/2 [J / 0 \ 2m+l 


ay 


K„ 


dy- 


s+l/2 


(|vWm) 


y 2 + cos 2 ip 


ay c ay . c 

ay sh — cos —h c ch — sin — 
2 2 2 2 


Then one may use the expansion 


= £■ 


sm 2k (f 


y 2 + cos 2 ip ( y 2 + l) fe+1 


which converges since y 2 + 1 ^ 1, sin 2 </? ^ l 13 . 


The next step is to integrate over y using the formulae (q ^ |r 

OO 

q 2u f dy I< v (q\fy 


2 i i \ / 2 /sh(n/)'l _ pr jrK v _ 3/2 {^ /g 2 - r 2 ) 
I ch (ry) j V 2 l K u _i /2 (\/q 2 - r 2 ) 


(C.l) 


(C.2) 


13 Notice, that the series is convergent even for the limiting value cos (p = 0, while the integral (|C.1|) converges for y —> 0, as a 
consequence of 

sh(au/2) a csin(|c|/2) a 2 

lim — — - 7 -— = - , lim —— 17 7 = — . 

y—>0 y 2 c—>-0 y 2 + c 2 /a 2 2 
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together with the identity 



(C.3) 


Explicitly, applying (IC.3I) once, the leading term is the one containing sh(ay/2) and is given by: 



OO 


2 m+2 r(m+ 1) 


c 

cos - . 
2 


0 


The first subleading term is 0(a~ 3 ): it comes from ch(ay/2) times basic Macdonald; the next subleading 
terms are 0(a -4 ): they come (i) from y sh(ay/2) times previous shift-index term; and (ii) ysh(ay/2) times first 
correction in sin 2 ^j. The end result is: 



[1] D. Gal’tsov, P. Spirin and T. N.Tomaras, JHEP01 (2013) 087, arXiv: 1210.6976 [hep-th]. 

[2] D. V. Gal’tsov, G.Kofinas, P. Spirin and T. N.Tomaras, Phys. Lett. B683 (2010) 331, arXiv:0908.0675 [hep-ph]. 

[3] V. Ferrari, P. Pendenza and G.Veneziano, Gen. Rel. Grav. 20, 1185 (1988). 

[4] G. Dvali, C. Gomez, R. S. Isermann, D.Liist and S. Stieberger, Nucl. Phys. B893 (2015) 187-235, arXiv: 1409.7405 
[hep-th]. 

[5] S. Weinberg, Phys. Rev. 140 (1965) B516. 

[6] S. Weinberg, “Gravitation and Cosmology”, Chapter 10; J. Wiley and Sons, 1972. 

[7] P.D.D’Eath and P.N. Payne, Phys. Rev.' D46, 658, 675, 694 (1992). 

[8] C. Herdeiro, M. O. P. Sampaio and C. Rebelo, JHEP 1107 (2011) 121, arXiv: 1105.2298 [hep-th], 

[9] Y. Const ant inou and A. Taliotis, JHEP 1311 (2013) 175, arXiv: 1308.2544 [hep-th]. 

[10] A. Gruzinov and G.Veneziano, arXiv: 1409.4555vl [gr-qc]. 

[11] T. Dray and G.’tHooft, Nucl. Phys. B253 (1985) 173. 

[12] R. Glauber, “Quantum Theory of Optical Coherence”, chapter 2; J. Wiley, 2007. 

[13] I. M. Gel’fand and G.E. Shilov, “Generalized Functions: Properties and operations”, Academic Press, Waltham, 
MA, 1964. 

[14] L.Smarr, Phys. Rev. D15, 2069 (1977). 

[15] P.C. Peters, Phys. Rev. Dl, 1559 (1970). 

[16] S. J. Kovacs and K.S. Thorne, Astrophys. J. 224, 62 (1978). 










